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We exploit exact inequalities that reler to the entropy of a distribution to derive a simple vari- 
ational principle at finite temperature for trial density matrices of Gutzwiller and Jastrow type. 
We use the result to extend at finite temperature the Gutzwiller approximation, which we apply to 
study a two-orbital model that we believe captures some essential features of V2O3. We indeed find 
that the phase diagram of the model bears many similarities to that of real vanadium sesquioxide. 
In addition, we show that in a Bethe lattice, where the finite temperature Gutzwiller approximation 
provides a rigorous upper bound of the actual free energy, the results compare well with the exact 
phase diagram obtained by the dynamical mean field theory. 
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I. INTRODUCTION 

A genuine Mott insulator, where the insulating char- 
acter is due exclusively to charge localization, is a very 
useful concept, physically conceivable, but never real- 
ized in the ground state of known correlated materials. 
Indeed, no system can sustain at zero temperature the 
residual entropy that would be associated with all other 
electronic degrees of freedom different from charge. As 
a result, Mott localization is always accompanied at low 
temperature by other phenomena that freeze those de- 
grees of freedom, for instance magnetic ordering or struc- 
tural distortions, which effectively turn the Mott insula- 
tor into a conventional band insulator. By this we mean 
the possibility of reproducing low-temperature static and 
often also dynamic properties of a supposed Mott insula- 
tor by an independent-particle scheme, no matter how 
sophisticated it isr^ However, even though it provides 
satisfactory results, an independent-particle scheme, like 
Hartree-Fock or DFT within LDA and its extensions, has 
a drawback: it can describe only the simultaneous lock- 
ing of charge and other degrees of freedom, like spin or 
lattice, while in a Mott insulator the charge freezes at a 
much higher energy scale than any other degree of free- 
dom. 

A tool that can reveal this hierarchy of energy scales 
typical of a Mott insulator is the temperature, which un- 
veils the profound difference between the excited states 
of a Mott insulator and those of its "band-insulator" 
counterpart^. 

The typical example are antiferromagnctic Mott in- 
sulators, e.g. La2Cu04 or Cr-doped V2O3, whose con- 
ducting properties remain insulating-like also above the 
Neel temperature, while any independent particle scheme 
would predict metallic behavior as soon as Neel order 
melts. It is therefore important to have at disposal 
theoretical tools able to deal with strong correlations 
and finite temperature. One that is currently adopted 
is dynamical mean field theory (DMFT)^ in combina- 



tion with local-density approximation (LDA)^ or the GW 
approximation^,. These combined methods are extremely 
reliable for correlated materials, but they can become ex- 
tremely cumbersome and numerically demanding, espe- 
cially if full consistency on the density is required. 



More recently, there have been several attempts to 
combine efficiently and self-consistently LDA with the 
Gutzwiller variational approach within the Gutzwiller 
approximation (GA),— which is less accurate but much 
simpler and less demanding than DMFT. So far, these at- 
tempts have been restricted to ground state properties, 
hence to zero temperature. A finite temperature exten- 
sion of the GA has actually been proposed in Ref. 0, 
but only in the simplest case of the one-band Hubbard 
model, where it was shown to reproduce qualitatively 
well the exact DMFT phase diagram<^ However, in or- 
der to tackle more realistic models, it would be desirable 
to have at disposal a finite temperature generalization of 
the GA able to deal with generic multi-band models. 



This is actually the aim of the present work that is or- 
ganized as follows. First, in section|lT]we derive a rigorous 
upper-bound estimate of the free-energy of a many-body 
Hamiltonian within the class of Gutzwiller- and Jastrow- 
like variational density matrices. Next, in section Hill we 
specialize to the case of Gutzwiller-like density matri- 
ces and introduce the Gutzwiller approximation at finite 
temperature. In section ITVl we introduce a simple model 
for vanadium sesquioxide, V2O3, show how to implement 
in practice the Gutzwiller approximation at finite tem- 
perature, section |IVA[ and discuss its results at zero and 
finite temperature, , section ITV Bl and flV C[ respectively. 
In section |V] we compare the phase diagram obtained by 
the Gutzwiller approximation on a Bethe lattice with the 
exact one obtained by DMFT. Finally, section IVll is de- 
voted to concluding remarks. 
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II. VARIATIONAL ESTIMATE OF THE FREE 
ENERGY 

In this section we shall repeatedly use some known 
trace inequalities, for which we refer the reader to Ref. 
[Tol . Let us consider an interacting many-body system de- 
scribed by the interacting Hamiltonian % at finite tem- 
perature T > 0. It is known that the free-energy func- 
tional 



F{X) = Tr(X-H) +TTr(XlnX), 



(1) 



with the matrix X > Q and such that TrX = 1, is mini- 
mized by the Boltzmann distribution function 



Trc'^«' 



(2) 



where /3 = 1/T. Therefore, any variational ansatz for the 
density matrix X provides an upper bound of the actual 
free energy 

F EE F(X,nin) < > with Tr X = 1 . (3) 

It is also known that, for any positive matrix y, the 
entropy of the distribution X satisfies the inequalityiS 

S{X) = -Tr(XlnX) > -Tr(Xlny) 

-i:i(x\n{XY-^))=S,,,{X,Y). (4) 



It then follows that, for any positive Y and X such that 
TrX = 1, 

F < min{Tr(XH) - T 5var(^, y)}- (5) 

Eq. ([5|) provides a variational principle for the free energy 
in terms of the distribution X and the matrix Y > 0. Let 
us assume the variational ansatz 



X = Vp,V\ 



where 



P* = 



Tre 



(6) 



(7) 



is the Boltzmann distribution corresponding to a varia- 
tional non-interacting Hamiltonian "H*, and V a many- 
body operator that we can parametrize as 



(8) 



with unitary U and Q > 0. It follows that the entropy of 
the distribution X 



Six) = -Tr(^XlnXj 

= -Tr(Qi/2p,Qi/2in(Qi/2^^gi/2)j^ (g) 



is independent of the unitary operator U. By means of 
Eq. dl]), setting F = Q, we obtain 



S^^AX, Y) = -Tr(^Qi/2 Qi/2 

-Tr(Qi/2p, Qi/2 ln(Qi/2p^ g-i/2) 
= -Trfp, Q In q) - Trfp, Q In (p,)). (10) 



In conclusion, given the ansatz Eqs. dn])-®, one can 
obtain an upper estimate of the actual free energy 

F < mmlTi-(^p^V'^nr^ +TTi-(^p^V''V InP^p) 

-TTi(p,V^rinpA\, (11) 



minimizing with respect to a non-interacting Hamilto- 
nian "H* and a many-body operator V. This minimiza- 
tion is feasible only for particular choices of V. For in- 
stance, if P = 1, Eq. pT|) reduces to the well-known 
Hartree-Fock variational estimate of the free energy. An- 
other possibility is that T' is a two-body Jastrow factor, 
which can be handled by the variational Monte Carlo sta- 
tistical approachiii In the next section, we shall consider 
still another class of operators V, which can be dealt with 
analytically in the limit of infinite coordination lattices. 

We conclude by noting that, since Eq. pT|) is based 
on the lower boimd estimate Sva.T:{X,Y) of the entropy 
of the distribution X, Eq. (|10p . there is no guarantee 
that, within a particular class of Y , such estimate is pos- 
itive at any temperature, as the true entropy should be. 
Therefore, it is more appropriate to state that 

S{X) > M^x{s,,,iX,Y),Oy (12) 

We further mention that Eq. (jj]) is actually the p = 1 
case of the more general inequalityiii 

S{X) = -Tr(XlnX) > -Tr(Xlnr) 



--Tr(Xln(X^r-^)j, 



(13) 



which becomes an equivalence as p ^ 0. We cannot 
exclude that exploiting Eq. ([T3|) one could get a better 
but still manageable estimate of the entropy, though we 
did not succeed. 



III. THE GUTZWILLER APPROXIMATION AT 
FINITE T 



We assume a Hamiltonian 

^ = E E i^^J + H.c.) + E (14) 

i,j ab i 

defined on a lattice with coordination number z, and hop- 
ping parameters such that their contribution to the 
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total energy is well behaved also in the limit z — > oo. Hi 
includes on-site potential and interaction terms, while 
creates an electron at site i with quantum number a, 
which can include spin and orbital. Given a variational 
density matrix of the form as in Eqs. ^ and (O, i.e. 

PG = Vp,V^, (15) 

we take the operator P to be of the Gutzwiller type^^^ 
i.e. 



(16) 



where Vi acts on the Hilbert space at site i. Wc shall 
impose that 



(17) 
(18) 



where Ci is any single-particle operator at site i. The 
above conditions replace those at zero temperature^^"— 
and allow to analytically compute averages over the dis- 
tribution function pQ in the limit of infinite coordination 
number, z — > oo (the proof is exactly the same as in Ref. 
[igI . hence we shall not repeat it here). Specifically, when 
z oo, the two conditions ((T7)) and (US]) imply that the 
distribution po has unit trace, and that all the formu- 
las of Ref. [la hold with the only difference that, instead 
of averaging over a variational Slater determinant, one 
has to average over the variational non-interacting Boltz- 
mann distribution . Following Ref. [igI . we assume that 
identifies a local natural basis, with creation operators 
djjj, related by a unitary transformation to the original 
operators c|^, such that 



Tr(^p*dLdj^) Sapn°^, 



(19) 



where n^^ depends on the variational Hamiltonian 
and on the temperature. We also introduce local Fock 
states in the natural representation 



h;M> = n(4)"° |o), 

a 

whose local probability distribution is diagonal, 
Tr(^p* I i;{n}){i;{m} \ ) S{n}{m} P°{„} 



(20) 



We then parametrize the operator Vi as 

V. = '^m^} {P'{n}) M;r)(*;M|, (21) 

{n},r 

where ^i-^r{n} the components of a temperature- 
dependent variational matrix $i, and | i;r) are local 



basis states in the original representation. In terms of 
$i the Eqs. (dZl) and (HH]) read 



Tr 



Tr($l$^) 
did 



1, 



(22) 
(23) 



all other bilinear operators in (|23p having null average. 
If we discard for simplicity the possibility of supercon- 
ductivity, it follows that^ 

Tr(pG ^) = E E (CC dld,0 + H.c. 



where 



.aP _ nt ,ab p 

* ij ~ / y aa ''ij ^ 



ab 



with the renormalization factors 



,-1/2, 



T>f*I4*. d,^ 



(24) 



(25) 



(26) 



In other words, the average over pQ of the Hamiltonian 
is equal to the average over the uncorrelated distri- 
bution p* of a renormalized hopping Hamiltonian plus 
the sum of local terms that depend only on the varia- 
tional matrices $i. 

We next need to evaluate the entropy. We note that, 
in the z — >■ oo limit, and for any, even non-local, single- 
particle operator C 



Tri^p.V'^VCj = Ti- i^pX 
Since it also holds that 

Tr (p.-p^-p InV^V^ = E Tr (p^-p/P^ In^l^^ 



ETr 



$1$, In ( (Pj 



it follows that Eq. ([TU| reads, in the z ^ oo limit, 

^varlp*,*^*) = S{p,) 



5(p.)+E^H*.|l^° 



(27) 



where S'^<I>|(f>j ||P°^ is the relative entropy between the 

distribution ^l^.^ and the uncorrelated local distribution 
. In conclusion, the free energy can be upper estimated 
through 

F < minj E E {t% 4^,^ + H.c.)) 

^ t-J al3 
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i 

-rMax(5var(p*,$^$),0 



(28) 



hence one just needs to minimize the r.h.s. supplemented 
by the constraints ([22]) and ((23|). A possible route is to 
regard n^^ in Eqs. (|19p and (|23p as independent mini- 
mization parameters, and introduce two Lagrange multi- 
pliers terms 



i a/3 



Sap 



where the non-interacting potential V enforces Eq. ([T9]> . 
while /iiQ^ enforce Eq. (|23p . 

When S'var(p*, > 0, the saddle point with re- 

spect to the uncorrelated distribution p*, see Eq. (O, 
corresponds to identifying 

^*-EE(Cf,4rf,^+^-^-)+v, (29) 

ij' a/3 

SO that, once V is chosen so as to satisfy Eq. p^ . Eq. 
reads 



F < mini ^4^',^°] 



E 



Tr 



EE 

; min { F [<i> , n'> , , 



s 
''a/3 riia 



(30) 



where F^, is the free energy of non-interacting electrons 
described by the Hamiltonian in (|29p that depends 
on the variational matrices $/ and on the parameters 
n°Q, through the constraint (|19p and the Eqs. (|25p and 
(IMl). 

However, when T ^ 0, S'var(p*, 'I'^^) becomes very 
small and it may happen that the saddle point of Eq. (pS)) 
is only a relative minimum, while the actual minimum is 
obtained by being the projection onto the ground state 
of the Hamiltonian "H* in Eq. ((29)) . In this well 
as when Svs^r{p*,^^^) < 0, the variational estimate of 
the free energy coincides with that of the ground state 
energy, evidently a drawback of the entropy bound that 
we use. In our experience, this problem may arise only 
at very low temperature, in which case, although the free 
energy estimate is continuous, its derivative may not be 
so. This is evidently not a signal of a genuine first order 
transition but only a flaw of the method. However, since 
the temperature when this occurs is extremely low, the 
entropy contribution is nonetheless negligible, so is the 
error in the free energy. 



Minimization of F[^,n^ , /x] therefore provides an up- 
per bound to the actual free energy in lattices with in- 
finite coordination number z — > oo. Seemingly to what 
it is done at zero temperature, one can keep using the 
same free-energy functional also when the coordination 
number is finite, which can be regarded as the finite tem- 
perature extension of the Gutzwiller approximation.—^— 
We mention that, in the simple case of a one-band Hub- 
bard model, the free energy functional i^[$,n°,/L<] coin- 
cides with the expression derived by different arguments 
in Ref. i. 

In the next section we apply the outlined method to a 
simple two-band Hubbard model and test its results with 
exact ones obtained by DMFT. 



IV. A TOY-MODEL FOR V2O3 

The model we are going to analyze is inspired by the 
physics of V2O3. In this compound, the ions have 
two valence electrons occupying the conduction bands 
that originate mainly from the t2g atomic c?-orbitals. At 
high temperatures, V2O3 is a paramagnetic metal but, 
upon substituting V with Cr it can turn into a param- 
agnetic insulator The transition is first order and ends 
into a second-order critical point. At low temperature, 
V2O3 is instead an antiferromagnetic insulator. The Neel 
transition occurs at T/v ~ 170 K and is accompanied by 
a martensitic transformation from the high-temperature 
corundum structure to the low-temperature monoclinic 
oneJ^ As a result, the magnetic ordering is not a sim- 
ple G-type, as it could well be in a bipartite lattice, 
but, in the honeycomb-lattice basal plane, two bonds 
are antiferromagnetic and one is ferromagnetic There 
is wide consensus that the magnetic moment is formed 
by a spin S = 1— but it is also contributed by angular 
momentumf2i signaling a non-negligible spin-orbit cou- 
pling. Even though a reliable description of the antifer- 
romagnetic transition requires including electron-lattice 
and spin-orbit couplings, the main features of the phase 
diagram can be likely explained ignoring those additional 
complications. The trigonal field of the corundum struc- 
ture splits the t2g orbitals into a lower Cg doublet and a 
higher aig singlet. It is therefore tempting to conclude 
that the low-temperature insulator describes the two 
electrons in the doublet that, because of Hund's rules, 
are coupled into a spin-triplet. This conclusion is proba- 
bly not far from reality. Indeed, although the bare value 
of the crystal field splitting is too small in comparison 
with the bare conduction bandwidth^^ strong enough 
electronic correlations may reverse the situation and sta- 
bilize the insulating phascj^ This scenario has been ac- 
tually advocated to explain the phase diagram of V2O3 
on the basis of a DMFT-LDA calculation in Ref. [13, and 
seems supported by some experimental evidences.— In- 
deed, DMFT-LDA results have shown that the effective 
crystal field splitting Aoft between and aig orbitals is 
enhanced by correlations from its bare value A ~ 0.27 eV 
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to one four times larger, which increases as the strength 
of the electron repulsioni^ In addition, Aeff has been 
found to increase upon lowering temperature T, though 
only slightly?^ but, more importantly, it has been ob- 
served that the magnetic susceptibility of the increases 
substantially with lowering T, while that of the aig stays 
constants^l or even diminishes^ precursor signals of a 
magnetic instability that involves only orbitals. 

However, all calculations so far have not been pushed 
down to the Neel transition temperature to really uncover 
the proposed mechanism of a gradual depopulation of the 
aig-derived band and concomitant magnetic polarization 
of the e^-derived ones. Here, we would like to address 
this issue by exploiting the finite temperature technique 
described in section IhI on a simplified model that we be- 
lieve captures the essential physics. Instead of consider- 
ing three t2g orbitals split into two plus one and occu- 
pied on average by two electrons, we shall consider only 
two split orbitals occupied on average by one electron 
(quarter filling). In this way we miss the important role 
of Coulomb exchange, which forces the two electrons on 
the Cg doublet to lock into a spin triplet state and might 
bring about relevant incoherence effects^^ but the gross 
features of the phase diagram, in particular the interplay 
between temperature, crystal field splitting, correlations 
and magnetism, should be maintained. 

Specifically, we shall study the two-band Hamiltonian 
on a square lattice 



crystal field A^s caused by U, which eventually leads 
to antiferromagnetism once the highest band is emptied, 
can be described also within Hartree-Fock. Indeed, if we 
neglect magnetism and assume the variational mean-field 
ansatz 



5n, 



{n2i 



5n, 



then the Hartree-Fock energies of the orbitals are ei = 
-A f/ (3 - 2(5n)/2, and £2 = A J7 (3 2(5n)/2, so 
that the effective crystal-field splitting is, within mean- 
field, Aeff = A -I- [/ (5n > A. As U increases, Aoff grows 
hence the highest band depopulates until it becomes com- 
pletely empty. Beyond this point, only the lowest band 
remains occupied, specifically half-filled, which can lead 
to a Stoner-like antiferromagnetic instability, hence to 
an insulating state. In other words, as we mentioned in 
the Introduction, an independent particle picture, like 
Hartree-Fock, is indeed able to explain the occurrence 
of an antiferromagnetic insulating state at low tempera- 
ture. However, no matter how large U is, Hartree-Fock 
will predict this insulating phase to turn metallic above 
the Neel temperature T^r. On the contrary, we expect 
that, for T > T/v but U large enough, the phase should 
still be insulating, though paramagnetic. 



A. Finite-T Gutzwiller approximation at work 



a=l ka 



E 



- A (nii ~ n2t) + ^ + n2tY 



(31) 



where a = 1, 2 labels the two orbitals, ek = — 2i ( cos + 
cos ky^ is the standard nearest neighbor tight-binding en- 
ergy, U parametrizes the on-site repulsion and A > the 
crystal field splitting. We include an inter-orbital hop- 
ping 7i(; = —At' sin k,^ sin ky with a symmetry such that 
the local single-particle density matrix remains diagonal 
in the orbital indices 1 and 2, thus mimicking the ai^-e^ 
hybridization in the corundum phase of ¥2031^ We shall 
further assume a density corresponding to one electron 
per site. 

In spite of its simplicity, the model in Eq. ([3T|) re- 
produces qualitatively the actual behavior of V2O3. If 
A ~ t' ^ which we shall consider hereafter, the model 
describes a two-band metal for small U . However, for 
very large J7, we do expect a Mott insulating phase with 
the electrons localized mostly on the lowest orbital and 
antiferromagnetically ordered. Therefore a strong repul- 
sion U can turn the two band metal into a single-band an- 
tiferromagnetic insulator, the two-band analogue of what 
is predicted in V2O3. The question we would like to ad- 
dress here is the behavior at finite temperature. 

We first observe that the enhancement of the effective 



We can improve the Hartree-Fock description at finite 
temperature by the Gutzwiller variational approach of 
section |TT1 By our choice, even though the inter-orbital 
hybridization t' is finite, hence the two orbital can mu- 
tually exchange electrons, still the local density matrix 
is diagonal by symmetry. In other words, the natural 
basis, see Eq. (fTg]) . coincides with the original one; a 
great simplification in the calculations. Since we will 
search for simple two-sublattice Neel order, we can set, 
for any site i belonging to sublattice A, <I>i = <I>^ = <& 



and 71°^^ 



- n\^^ = (see Eq. HI])), such that 

1. The variational matrix $ is defined in 



the local Fock space and is only invariant under spin- 
rotations around the magnetization axis, which we choose 
as the z axis. It follows that, for any site i belong- 
ing to the other sublattice B, = ^ with 

U = exp {iTrSy/'l^ and Sy the y-component of the total 

local spin, while n°^^ = n^^^ = Because of Eq. 



we must impose the constraint 



Tr( p*c. 



^ c 

iaa iha' 



Tr($tci,4^c,^, 



(32) 





a — a' 



for i G A, while, for i G i3, $ — >■ J/'l' $ J7 and — >• n\ 

Since natural and original bases coincide, the renor- 
malization factors of Eq. (pS)) arc diagonal and read, for 
i e A, 



Tr 



.t , 

'aa 



(33) 
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while, for i € B, R*^^ = K-a- Wc find that, at the 
optimized values of the variational parameters, Ri aa a-re 
always real. Therefore, if we define Raa = Ra + o' Sa G 
5fte, then the variational uncorrelated Hamiltonian "H*, 
see Eq. is 



2.5 



a=l kiT 
C. 



lkcr^2kcr 



+ E '^^"k(4k^C2i^+Q<T + ^-C- 



k(T 



"2ko 



ako" ko" 

2 

- ^ ^ a /la (clk^Cak+Q^ + H.c.) , 
a— 1 kcr 



(34) 



where Q = (7r,7r), the Lagrange multipliers /i, /^cf, /ii 
and enforce the constraints ([5^ . and 

Eak = (i?a ~ '5'q) ek, 
7k = {Ri R2 + Si S2) iw, 

7"k = (i?1^2 + i?2^l)7k. 

It follows that, if 

F4$,n°] = -^lnTr(e-'3'^*), (35) 
where N in the number of sites, then we have to minimize 
F[A,$,ri°] 



^* + f Tr($t(„^ + 



-rs'($^$||p°) 



Tr( $^$1 



(36) 



with the constraint y^„„ = 1. 

We find more convenient to minimize the variational 
free energy p6p first with respect to all parameters except 
■nP^ thus obtaining the functional 



F[n"] = minF[A,$,n°]. 



(37) 



We calculate F[n'^] in a two-step cyclic process; first we 
fix <i> and minimize F*[$,n°] with respect to the La- 
grange multipliers in Eq. (j34p . Then, at fixed matrix 



elements (c,Lg.c w^,)p, , we minimize F[X,^,nP] with re- 
spect to $ fulfilling the Gutzwiller constraints. This sec- 
ond non-linear constrained minimization is performed by 
the LANCELOT B routine of the GALAHAD library.^ 
This two-step cycle is repeated until convergence. Fi- 
nally a full minimization of _F[n'^] with respect to nP is 
performed. 
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FIG. 1: (Color online) T = phase diagram for model 
in Eq. The solid line indicates the MIT transition; 

the dashed/dotted line separates the two-band paramagnetic 
metal (2-band PM) from the one-band metal (1-band PM). 
The vertical dashed red lines indicate the values of A which 
are used in Fig[2] 



B. T = phase diagram 

The results that follow are obtained setting t = 1/8 
and the inter-orbital hybridization t' = 0.3t. The t' = 
bandwidth VF = 8< = 1 hence sets the unit of energy. 

First we consider the T = case of Eq. ([36|) . which 
corresponds to the usual Gutzwiller variational approach. 
In Figlljwe plot the zero temperature phase diagram in 
the paramagnetic sector as a function of U and A. Our 
results compare well with the DMFT phase diagram of 
Refs. [2^ and [2^. In the limit A = the model under- 
goes a second order metal-to-insulator transition (MIT) 
at a critical value U^^^ ~ 2.27VF. In the opposite non 
interacting case, C/ = 0, upon increasing A the system 
crosses a Lifshitz transition from a two-band to a one- 
band metal. We note that the majority (>) and minority 
(<) bands do not have a unique orbital character, there- 
fore the band polarization n> — n<. which saturates to 
1 at the two-band — > one-band transition, is in general 
different from the orbital polarization ni — 77.2. 

As anticipated, at finite U a smaller crystal field spit- 
ting is required to induce the two-band — >■ one-band tran- 
sition, see FigHJ Above this transition, the ground state 
is a one-band metal which eventually undergoes a sec- 
ond order MIT at a critical value Uc — 1.68VF. In 
Fig. [2] (right panel) we show the details of these two 
subsequent transitions for a value of A = 0.3VF. For 
U < 0.64VF, the two-band metal is stable but, increasing 
[/, the minority band gradually empties and both renor- 
malization factors, i?i and R2, decrease. At ~ 0.64IF 
the minority band completely depopulates and the left- 
over half-filled majority band is driven to the MIT at 
U ~ 1.68IF. Approaching the MIT, the renormalization 
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A=0.025W A=0.3W 

^ I ' I p I I • I r \ ' \ r 
I I I 




FIG. 2: (Color online) n^—n^ (blue diamonds) , renormaliza- 
tion factors Ri (black squares) and R2 (red circles) as a func- 
tion of U for fixed A = 0.025W (left panel) and A = 0.3W 
(right panel). The vertical dashed lines indicate the two-band 
— >■ one-band metal transition and the MIT. 




2AAV 

FIG. 3: (Color online) Zero temperature phase diagram al- 
lowing for magnetism. The black line separates the two-band 
paramagnetic metal (2-band PM) from the one-band antifer- 
romagnetic insulator (1-band AFI). In the inset we plot the 
orbital polarization ni — 712 and the staggered magnetization 
for a fixed value of A = 0.0251^ (red dashed line). 



factor Ri of the lowest-energy orbital vanishes, while R2 
actually increases to one - the almost empty orbital un- 
dresses from correlations. For smaller A, the two-band — ?> 
one-band transition becomes first order and approaches 
the one-band MIT point, ending in a multicritical point 
at A ~ 0.0281^. In Figl2](lcft panel) we plot the behav- 
ior of Ri, i?2 and (n> — n<) for A = 0.025VF; in this case 
the two renormalization factors are approximately equal 
and decrease monotonically with U. At the transition, 
the majority orbital occupation suddenly increases and 
the corresponding renormalization factor vanishes. We 
mention that a Mott insulator with partial occupation of 
both orbitals can not be stabilized within the Gutzwiller 
approximation, while more reliable DMFT calculations 
show that such a phase does exist for very small 

If we allow for magnetism, the one-band phases, ei- 
ther metallic or Mott insulating, turn immediately into 
an antifcrromagnctic insulator, sec Figj^l The transition 
from the two-band paramagnetic metal (2-band PM) to 
the one-band antifcrromagnctic insulator (1-band AFI) is 
first order and accompanied by a jump in the orbital po- 
larization and in the staggered magnetization, see inset 
of FigEl 



C. T 7^ phase diagram 

We have seen that at zero temperature the ground 
state is either a one-band antiferromagnetic insulator or 
a two-band paramagnetic metal. Therefore, as we an- 
ticipated, the T = Gutzwiller variational results are 
not dissimilar from the predictions of the Hartree-Fock 
approximation. Differences instead arise at finite tem- 
perature, where the Gutzwiller variational approach, as 



we are going to show, can describe melting of the Neel 
order without metallization, unlike Hartree-Fock. 

Wc begin as before by restricting the analysis to 
the paramagnetic sector and consider the case of A = 
0.025W. At zero temperature we found that the model 
is a 1-band PI for values of J7 > 1.7M^, while a 2-band PM 
below, Fig[T] At finite temperature, the entropic contri- 
bution may favor the paramagnetic insulating solution, 
like in the single band Hubbard model;^ thus leading to 
a finite T metal-insulator transition. This indeed occurs, 
as shown in Fig. |4] where we plot the phase diagram as 
a function of U and T (upper panel) and the tempera- 
ture dependence of the majority orbital i?i and the or- 
bital polarization (lower panels). In the figure we observe 
that for values oi U > 1.7W, increasing the temperature 
the orbital polarization decreases and the quasiparticle 
weight increase: the 1-band PI continuously evolves to- 
wards a 2-band PI. Instead, for smaller values of U, the 
system is initially a two-band metal and undergoes a first 
order transition to an insulating state which is accompa- 
nied by an abrupt fall-down of the renormalization fac- 
tors and increase of orbital polarization. As in Ref. 
we interpret the jump of the renormalization factor as 
the boundary of the PM-PI transition. Notice that, dif- 
ferently from the T = case, the orbital polarization 
does not saturate at the transition. Finally, for values of 
U smaller than ~ 1.191^, the quasiparticle weight and 
the orbital polarization evolve smoothly to the high tem- 
perature limit, displaying a dip that we interpret as the 
crossover regime. We estimate the end-point of the tran- 
sition at T ~ 0.09W. 

Wc note, in the lower panel of Fig. U] and for U = 
I.IM^, the tiny discontinuity of Ri and ni — n2 at T ~ 
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O.OIM^, which is consequence of the aforementioned arti- 
ficial discontinuity in the slope of the free energy caused 
by our not rigorous lower bound of the entropy. 




TAV 



FIG. 4: (Color online) Upper panel: phase diagram in the 
paramagnetic domain. The black line separates the PM phase 
from the PI one. The red vertical lines indicate the values of 
U plotted in the lower panel. 

Lower panel: Temperature dependence of the quasiparticle 
renormalization factor for the majority orbital (left) and of 
the orbital polarization (right) for different values of U. 

When magnetism is allowed, at zero temperature and 
at large U the ground state is antiferromagnetic. At finite 
temperature the system remains ordered up to the Ncel 
temperature. In Fig. [5] we plot for A ~ 0.0251^ the 
phase diagram, indicating by a dotted line the PM-PI 
transition that we have found in the paramagnetic sector. 
We note that the PM-PI transition line crosses the Neel 
temperature, roughly at C/ ~ 1.281^, and extends above. 
For U > 1.281^, the Gutzwiller variational approach is 
able to describe melting of the AFI into a 2-band PI, 
which we mentioned is not accessible by Hartree-Fock. 




UAV 

FIG. 5: (Color online) Finite temperature phase diagram as 
a function of U and T, for a fixed value of A = 0.02511^. 
The paramagnetic solution is continued also within the AFM 
domain (dotted line). The vertical red dashed line indicates 
the temperature cut represented in Fig[6l 



For smaller values of C7, the magnetic insulator turns into 
a 2-band PM that eventually undergoes a Mott transition 
at higher temperatures. In Fig. IHlwc show more in detail 
the behavior of the physical quantities across the different 
transitions; in the low temperature AFI (blue area on 
the left) the orbital and magnetic polarizations are very 
weakly temperature dependent. Increasing T, the Neel 
order melts, the orbital polarization decreases (red areas 
on the right), and the model turns into a 2-band PM 
(left panel) or PI (right panel) depending on the value of 
U. In the former case, left panel of Fig. [Bl the 2-band 
PM is eventually driven to the PI state, transition that is 
signaled by the sudden vanishing of the renormalization 
factors and the jump of the orbital polarization ni — 712. 
In the right panel, instead, the AFI melts directly in the 
PI; the renormalization factor vanish at the transition 
and then smoothly increases from zero on raising T. 

We observe that the finite temperature phase diagram 
of Fig. [5] is not dissimilar to that of V2O3 as function 
of chemical/physical pressure. Also the physical mech- 
anism that controls the phase diagram, i.e. the corre- 
lation enhanced crystal-field splitting, is consistent with 
that proposed in Ref. [13 for V2O3, though in our case 
the number of orbitals involved is two and not three. We 
also note the discontinuous increase of the orbital polar- 
ization across the PM to PI transition upon increasing 
temperature, see right panel in Fig. [6l which is consis- 
tent with X-ray adsorption spectra of V203 )^°'^° in which 
case the orbital polarization relates to the occupation of 
the Cg orbitals with respect to the aig one. 
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FIG. 6: (Color online) The blue (to the left in each panel) and 
red areas (to the right) indicate respectively the AFI phase 
and the paramagnetic phases as a function of temperature at 
fixed values of [/ = 1.22W (left panel) and U = 1.8W (right 
panel). At low temperatures the orbital polarization (blue 
squares) and the staggered magnetization (black squares) are 
practically equal to the zero temperature values and display 
a discontinuous jump at the AFI-paramagnetic transition. In 
the paramagnetic phase we show also the behavior of the 
renormalization factors whose jump indicates the PM-PI tran- 
sition. 



V. COMPARISON WITH DMFT 

In this section, we compare the quality of the fi- 
nite temperature Gutzwiller approximation with exact 
DMFT results. In particular, we shall consider a sim- 
plified version of the model in Eq. ()31|) with vanishing 
inter-orbital hybridization, t' ~ 0, and on a Bcthe lat- 
tice with only nearest neighbor hopping, which leads to 
a non-interacting semicircular density of states of band- 
width W. We choose a Bethe lattice (a Cayley tree with 
coordination number z — oo) because in this case DMFT 
is exact and, as previously discussed, the Gutzwiller ap- 
proximation does provide a rigorous upper bound to the 
free energy, which therefore makes it possible to assess 
its accuracy with respect to exact results. The phase 
diagrams obtained by DMFT and by the Gutzwiller ap- 
proximation in the U-T space for A = 0.025W are shown 
in Fig. [71 

DMFT maps the lattice model onto an impurity model, 
which, in the present calculation, is solved by means of 
exact diagonalization^ in the finitc-tcmpcrature imple- 
mentation proposed in Ref. [3^ which is particularly ac- 
curate at the low temperatures that we consider. Within 
the exact diagonalization approach, the bath is approxi- 
mated by a finite number, Nh, of energy levels. Here we 
take Nh = 10 and iVf, = 12 , i.e. 5 and 6 bath levels 
for each physical orbital. Only for Ni, = 10 we could 
include a number of states sufficient to obtained con- 
verged results. Therefore data for Ni, — 12 have only 
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FIG. 7: Top panel: finite T phase diagram within the 
Gutzwiller approximation at A = 0.025Vy, t' — and a 
semicircular DOS. Bottom panel: same as before but within 
DMFT, which is exact. 



obtained for low temperatures, and used to prove that 
the discretization error only leads to minor corrections 
to the phase diagram. We consider both paramagnetic 
and antiferromagnetic solutions. As customary, we first 
determine the Mott transition line in the paramagnetic 
sector by comparing the free energies of the metallic and 
Mott insulating solutions. The transition is first-order at 
any finite temperature and ends in a finite-temperature 
critical point at T < 0.05VF. If we allow for long-range 
antiferromagnetic order, at T = the system is metallic 
for U '^W , and it turns into a single-band antiferromag- 
net for larger values of the interaction. The Neel tem- 
perature rapidly grows with U and reaches a maximum 
around U ~ 1.2W, above which it monotonically de- 
creases. However, differently from the single-band case, 
the first-order Mott line is not completely covered by the 
antiferromagnetic dome. 

The DMFT phase diagram is thus very similar to that 
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obtained by the finite-temperature Gutzwiller approach, 
quahtatively and to same extent also quantitatively, see 
Fig. [71 As common with the Gutzwiller approxima- 
tion, the T = Mott transition in the paramagnetic sec- 
tor occurs at larger U/W ~ 1.7 than the exact DMFT 
value U/W ~ 1.5. In addition, the Gutzwiller wavefunc- 
tion seems to overestimate antiferromagnetism, which 
occupies a larger region in the phase diagram. How- 
ever, quite remarkably, the critical endpoints of the PM- 
PI Mott transition do not differ much, U/W ~ 1.17 
and T/W ~ 0.07 in the Gutzwiller calculation, while 
U/W ~ 1.15 and T/W ~ 0.05 in DMFT. 

VI. CONCLUSIONS 

Using some rigorous trace inequalities, we have de- 
rived an upper-bound estimate of the free energy of an 
interacting-electron Hamiltonian for variational density 
matrices of Gutzwiller and Jastrow type. We have then 
exploited this result to extend to finite temperature the 
conventional Gutzwiller approximation, which in turn be- 
comes an exact variational approach in lattices with in- 
finite coordination number. 

We have applied this technique to calculate the finite- 
temperature phase diagram of a two-band model that we 
believe captures qualitatively well the physics of vana- 



dium sesquioxide, V2O3. In spite of being extremely sim- 
plified with respect to a complete description of V2O3, 
the model has a very similar phase diagram comprising a 
low-temperature antiferromagnetic insulating dome and 
high-temperature paramagnetic metal as well as Mott 
insulating phases separated by a first order line with 
a second-order critical endpoint. We have tested the 
accuracy of our finite temperature Gutzwiller approxi- 
mation comparing the phase diagram of the model on 
a Bethe lattice with the exact one obtained by DMFT. 
The agreement is qualitatively very satisfying and partly 
also quantitatively. We believe therefore that this simple 
variational technique is very promising to attach corre- 
lated electron systems at finite temperature, and could be 
used whenever more reliable tools, like DMFT, become 
numerically too demanding. 
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